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Abstract. This is a report on recent results with A. Hassell on quantum 
ergodicity of boundary traces of eigenfunctions on domains with ergodic bil- 
Uards, and of work in progress with Hassell and Sogge on norms of boundary 
traces. Related work by Burq, Grieser and Smith-Sogge is also discussed. 

1. Introduction 

This is a report on recent work (and work in progress) with A. Hassell and 
C. Sogge on boundary traces of eigenfunctions of the Laplacian on bounded do- 
mains and their relations to the dynamics of the billiard map. Boundary traces of 
eigenfunctions arise when one makes a reduction of the interior/exterior eigenvalue 
problem or wave equation of a domain to its boundary. In the classical approach 
to boundary problems of Neumann and Fredholm, the boundary reduction is made 
via layer potentials and associated boundary integral operators. Our results are 
based on an analysis of these operators as semiclassical Fourier integral operators. 

We work on a Lipschitz domain C M" which is assumed to be compact and 
piecewise smooth with corners. We give a precise definition below. The interior 
eigenvalue problem is: 

-Auj ^ Xjuj in Q, {uj,Uk)=6jk 



Buj = on dfl, 



■i-jldn (Dirichlct), resp. 



(1.1) 
where 

(1.2) Bu, = 

dyUj + A'ujgn, (Neumann, Robin) 

Here, dv denotes the unit interior normal derivative. 

By the boundary trace it^ of an interior eigenfunction Uj, we mean the comple- 
mentary operator to _B, namely 

(1.3) u^jiq) — dyUj\dQ, (Dirichlet), resp. u^j{q) — Uj\dQ, (Neumann, Robin). 

These boundary traces are a reduction to the boundary of the eigenfuncions, and 
they become eigenfunctions of a reduction to the boundary of the resolvent, given 
by the boundary integral operator F\ on dU, defined by 



-2 /an ar7Go(y, 2/', X)f{y')dA{y'), Neumann 



(1.4) F^f{y) = 



'2lon^Go{y,y\X)f{y')dA{y'), 



Dirichlet 
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where dA{y) is the induced boundary area form, and where 

Go{z,z',X) = (2^)-"A"-2/e^M---')-«^,^de 

= CA"-2(A|z - z'|)-("-2)/2 Ha|^/^2_^(A|z - z'\). 

is the free outgoing Green's function K". Here, HeS^j^^iJ the Hankel function. 
Thus, the Dirichlet boundary integral operator has kernel 

where Vz is the unit inward pointing normal at z. It is elementary and classical 
(from Green's formula) that 

{Up Neumann 
— u^, Dirichlet. 

One may think of -Fa as a transfer operator or effective Hamiltonian in the reduction 
to the boundary. More precisely, as we will explain in Section |31 Fh is the quanti- 
zation of the billiard map /3 on the ball bundle B*dft, and is thus the reduction to 
the boundary of the quantum dynamics of the wave group in the interior. Yet it is 
clearly a more elementary operator than the Dirichlet and Neumann wave groups. 

Our interest is in the relation between the dynamics of f3 on B*dft and the 
asymptotics of the boundary traces -uj. In particular we study: 

• Asymptotics of matrix elements {A{hj)u^, u^); 

• Asymptotics of i'' norms: HmJUlp. 

• Asymptotics of ratios | |i,p/| |uj| 1^2 

The problems and results we present are to some extent analogues for boundary 
traces of results of |SZ[ IZZwj on eigenfunctions on boundaryless manifolds or on 
interior eigenfunctions of manifolds with boundary. In the next section, we state 
the main results at this time. 

Besides A. Hassell and C. Sogge, the author would like to thank M. Zworski for 
helpful discussions about non-convex domains and ghost orbits, and D. Tataru for 
checking the statements of his results in Section O (they are stated differently in 

ID) 

2. Statement of results 

Let be a bounded subdomain of R" with closure fi. 

Definition: We say that il C M" is a piecewise smooth manifold if the boundary 
Y = dn is strongly Lipschitz, and can be written as a disjoint union 

where each Hi is an open, relatively compact subset of a smooth embedded hyper- 
surface Si, and S is a closed set of (n — l)-measure zero. 

The sets Hi are called boundary hypersurfaces of fl. We call S the singular set, 
and write Y° = Y\ Y^ for the regular part of the boundary. 
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2.1. Statement of results on quantum ergodicity. Our main results concern 
the quantum ergodicity of boundary traces. To state the results, we need some 
notation. In the table, k denotes a C°° function on Y while k is the principal symbol 
of the operator K e '^^{Y), and da is the natural symplectic volume measure on 
B*Y. We also define the function 7(17) on B*Y by 

(2.1) l{q) = ^Jl-\ri\\ q={y,v)- 



Boundary Values 


B 


Bu 


u" 


dfiB 


Dirichlet 


U Y 


dyU Y 


j{q)da 


Neumann 


dyU\Y 


My 


lilT'da 


Robin 


{dyU — Ku) y 


My 


lilT'da 


-Robin 


{d^u - Ku)\y 


My 


l{q)da 





Theorem 1. Letfl C M" be a bounded piecewise smooth manifold (see Definition\^ 
with ergodic billiard map. Let {u^} be the boundary traces f the eigenf unctions {uj} 
of Ab on L^(f2) in the sense of the table above. Let Afi be a semiclassical operator 
of order zero on Y . Then there is a subset S of the positive integers, of density 
one, such that 

lim (Aii.u^,,u'') —LUBiA), B= Neumann, Robin or '^^ -Robin, 

(2.2) 

lim \-'^{Ahu],u]) = lob{.A), Dirichlet, 

where hj = XJ^ and ujBiA) = ^^.^^aAd^B- 

We note that the boundary traces are only normalized by HMjULaj-Q-) = 1. We 
obtain new results on L^ norms by taking A = I . For the Neumann case, we have 



(2.3) lim ||u 



, 2 _ 2vol(y) 



while for the Dirichlet boundary condition they imply 
(2.4) lim AT^IIu 



2 fc||2 _ 2vol(y) 



To our knowledge, the only prior results on L^ norms of boundary traces are 
the upper and lower bounds of jBLRI IHT| in the Dirichlet case, and the bounds 
implicit in T^ in the Neumann case. We review the latter in the last section. 

2.2. Statement of results on norms of boundary traces. We now consider 
L^ norms of boundary traces for general domains. We state a number of results of 
work in progress |HSZ| . For simplicity of exposition, we only state the results for 
norms. 

First we consider estimates of boundary traces and state general L°° bounds anal- 
ogous to those which follow from the remainder estimate of Avakumovic-Levitan 
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for the spectral function of the Laplacian on a boundaryless manifold (see | SZ| for 
references). We denote by 

(2.5) eIm^x,x)^ 

j:\j<Xe[a,b] 

the boundary trace of spectral projections kernel (on the diagonal) for VA corre- 
sponding to the interval [a, 6]. 

Proposition 2.1. We have: 

j:\j<X 

Hence, 

( 0(A"+i), Dirichlet 

II4a,a+i](^'^)IU~ = { 

[ 0(A""i), Neumann. 

Of course, the estimate holds for each term ||u^(9)P||L~(aa) separately. We state 
it in the above form because it is sharper than the statement for boundary traces 
of eigenfunctions: for instance, the results implies the same bound for boundary 
traces tpj of quasimodes of order satisfying ||(A — A|)'!/'j|| = 0(1). We also wish 
to emphasize that our estimates on boundary traces of eigenfunctions cannot be 
better than estimates for the above spectral projections. 

Even for boundary traces of eigefunctions, the above estimate on L°° norms 
is sharp among all compact Riemannian domains with boundary. Indeed, they 
are achieved on the northern hemisphere of the standard S" by zonal spherical 
harmonics with pole on the boundary: let xq be a point on the boundary of the 
northern hemisphere (i.e. the equator), and let Yf^^° be the L^-normalizcd zonal 
spherical harmonic of degree N with pole at xq, i.e. the spherical harmonic which 
is invariant under rotations fixing xq. Then ^[Y^^''{x) +Yf^^°(x*)] is a Neumann 
eigenfunction which is essentially i^-normalized and which has the same norm 
as Yf^^°. Here, x* is the reflection of x through the equation, i.e. (xi, . . . ,x„)* — 
(xi, . . . , Xn-i, —Xn)- Taking the normal derivative of the odd part of Y^ under the 
reflection produces a similar sharp example in the Dirichlet case. Eucliean examples 
will be discussed in the final section. 

However, our main result on norms of boundary traces, as in the interior case 
studied in ISZI , shows that this result is generically not sharp. In the following, 
the symbol f{x) = fl{g{x)) is the negative of f{x) = o{g{x)), i.e. f{x) grows at 
least as fast as g{x) along a subsequence. In the following, -E[a,a+o(i)] is short for 
E'[a,b] where a, b are the endpoints of a sequence of shrinking intervals satisfying 
a = A, = A + e(A) where < e(A) = o(l). 

Theorem 2. Suppose that 

{^iXj^), Dirichlet 
57(Aj ^ ), Neumann. 



CA"+2 + 0(A"+i), Dirichlet 
A" + 0(A"-i), Neumann. 
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or more generally that 

( Dirichlet 
\\E[x,x+o{i)\{x,x)\\l^ = I 

y fi(A~2— ). Neumann. 

Then there exists a recurrent point q G dfl, i.e. there exists T > and a positive 
measure set of directions e S*fl such that the billiard orbit exp^(i^) returns to q 
at time T. 

To put this another way, if the set of looping directions at q has measure zero 
for all q S dfl, then for any e > 0, 

(2.6) limsupA"'T^||£;[Aa+e](9,'?)l|L- ^0 as e 0. . 

A — ^oo 

As pointed out in Section [5.4l a recurrent point in the sense above cannot occur in 
a convex analytic Euclidean domain. Hence, 

Corollary 2.2. // C M" is a convex analytic domain with Euclidean metric, 
then 

{ "+1 
o(A^ ^ ), Dirichlet 
o{Xj ^ ), Neumann. 

The hypothesis in the Corollary can probably be weakened quite a bit. For 
instance, it seems unlikely that real analytic Euclidean domains of any kind can 
have recurrent points in our sense, and the same may be true for smooth Euclidean 
domains (see also fSV, for related conjectures). 

We note that although the result on boundary traces appears similar to that of 
|SZj on eigenfunctions on manifolds without boundary, there is a serious difference 
in that the estimate in the boundaryless case is on the ratio | |wj | |lp/| |uj | II^ whereas 
in the case of boundary traces, it is the interior eigenfunction and not the boundary 
trace which is normalized. In the Dirichlet case, | |l2 is known to be bounded 
above and below by a constant times Xj BLR, HT', so one could restate the result 
in terms of the ratio ||u^||ioo/||w^||^2. In the Neumann case, ||mJ||^2 can vary with 
the domain and metric. 

2.3. Remarks on methods and related results. As mentioned above, our ap- 
proach to boundary problems is to work entirely on the boundary. In a recent article 
jBuj , Burq takes the opposite route of 'reducing' the study of ergodicity of bounday 
traces to the interior. He extends the methods of Gerard-Leichtnam |GH to reduce 
the proof of ergodicity of boundary traces to the results of Zelditch-Zworski |ZZw| 
on ergodicity of interior eigenfunctions. Both approaches have their advantages, to 
which we devote a few remarks. 

Our motivation for working entirely on the boundary comes from several sources. 
First, it seems to us natural to prove results about boundary traces by working with 
quantum dynamics on the boundary, i.e. by making a semiclassical analysis of the 
operator Fh- We also hope (and believe) that the semiclassical analysis of Fh and 
related operators has an independent interest and other applications, for instance 
to inverse spectral theory (see jZlp. 

A second motivation is that the use of semiclassical asymptotics of the resolvent 
and of Fh to relate spectrum and billiards is now the standard approach in the 
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physics literature on eigenvalues of bounded domains. The method originated in 
the articles of Balian-Bloch BBl BB2 and has been developed in, for instance, the 

articles |XnilB> .BS. .BFSS GP . THS..TH;?21ItV| . A key reason for the wide use of 

the boundary reduction (as explained to the author by A. Backer and R. Schubert) 
is that boundary traces of eigenfunctions on smooth domains are much easier to 
compute numerically than interior eigenfunctions. See the expository article 
by A. Backer (especially sections (3.3.1) - (3.3.7) and figure 20) for a discussion of 
the numerical methods, for further references, and for pictures of boundary traces 
of eigenfunctions. Despite the classical nature of layer potentials and the induced 
boundary integral operators, and their common use in physics and engineering, we 
are not aware of prior mathematical studies of their semiclassical properties. 

The boundary integral method based on seems to us more elementary than the 
interior methods, but it does have the disadvantage that the phase of Fh, given by 
the distance function between boundary points, is the generating function of both 
the interior and exterior billiard maps. Hence one has to deal with the so-called 
ghost billiard orbits which have links outside the domain (see Section l4.1.4|l . In the 
original version of JHZJ , we restricted to convex domains where ghost orbits do not 
arise. Using his interior method, Burq B then proved a more general ergodicity 
result for non-convex as well as convex domains. After that, in the second version 
of |HZ| . we extended our methods and results to general (possibly non-convex) 
domains. In the interim, discussions with M. Zworski were very helpful. The 
additional complication of ghost orbits of non-convex domains does not in the end 
turn out to be serious. 

3. Boundary integral operators as semiclassical Fourier integral 

operators 

We begin by explaining the sense in which Fh is the boundary reduction or 
effective Hamiltonian for the resolvent. It is a classical result of potential theory. 

We recall that the single, resp. double layer potentials of a domain Vi C are 
the operators 

r Sl{\)f{x) ^ jg^G^{x,q,\)f{q)dA{q), 

(3.1) <^ 

{ ViiX)f{x) = /g,, g§-GoX, q, X)fiq)dAiq), 

where dA{q) is the surface area measure on dfl, where v is the interior unit normal 
to n, and where 9^ = ;/ • V. They induce the boundary integral operators H1.4() . 

Denote by i?n(A) (resp. i?Qc(A)) the resolvent of the interior Dirichlet Laplacian 
(resp. the exterior Neumann Laplacian). Then the Fredholm-Neumann reduction 
takes the form: 

r Rf,{X) = loi?o(A)ln + InmMl + Fxr'SiXy^n, 

(3.2) <^ 

[ R^.'^iX) = loe<(A)loe + ln^VeiX){I + Fx)-'SiiXy-ln^ 

This formula has the form of a Grushin reduction of the Laplacians to operators 
on the boundary in the sense of Sjostrand (see e.g. |DSj ). Although we are not 
studying trace formula here, we note that the combination of these two formula 
implies that (in the sense of distributions), 

TrM2[R%\X) (BR%~ Rop{X)] = ^ logdet(/ + F^), . 
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where det(/ + F\) is the Fredholm determinant. This formula is often used to 
determine the eigenvalues and resonances of domains. 

It should be observed that the Dirichlet resolvent has the complexity of the 
Dirichlet wave group, while Fx is no more complicated that the free resolvent. It 
is {I + F\)~^ which has the complexity of the Dirichlet resolvent or the Dirichlet- 
to-Neumann operator. The relative simplicity of F\ is one of the attractions of the 
boundary reduction. 



3.1. Quantized billiard map: Neumann. Since Fx is a reduction to the bound- 
ary of the resolvent, it should be a kind of quantization of the billiard map. We 
now explain in what sense this is true. 

The billiard map /3 : B*F° f*Y is defined on the open ball bundle B*Y° as 
follows: given (y, 77) S T*F, with |?7| < 1 we let (y, C) S S*U, be the unique inward- 
pointing unit covector at y which projects to (y,?]) under the map TqqP^ — > T*Y . 
Then we follow the geodesic (straight line) determined by (y, to the first place 
it intersects the boundary again; let y' €Y denote this first intersection. If y' £ E 
then we define (3{y, r]) — y' . Otherwise, let 77' be the projection of C, to T*,Y . Then 
we define 

The map /3_ : B*Y° — > T*Y is defined similarly, following the backward billiard 
trajectory (that is, the straight line with initial condition (y, 2(C • i'y)vy — C))- 
The graph 

(3.3) Cbiu = graph/3 = {{I3{q), q)\qe n"}. 

of /3 is a smooth Lagrangian submanifold of B*Y° x B*Y° . In a neighbourhood of 
(yo, Vo,y'o,v'o) it is given by 

(3.4) Cbiu = {{y, -Vyd{y, y'), y' , Vy.d{y, y')), 

where d{y,y') is the Euclidean distance function. For strictly convex it is given 
globally by (|3.4|) . for y,y' G Y°, but this is not true in general. This causes extra 
difficulties for nonconvex domains, namely it introduces spurious billiard orbits 
(known as ghost orbits in the physics literature) which do not remain entirely in 
the domain, and which have to be shown to be irrelevant. 

We see directly from (|1.4ll that Fh is an oscillatory integral operator with phase 
equal to the generating function d{y, y') of Cbiu- This is especially clear in dimension 
3 when 

-j^ ^i\\x — x'\ 

G[x, X ,\) = - — j — . 

47r \x — x'\ 

However, it is also clear that Fh{y,y') has the homogeneous singularity on the di- 
agonal of a pseudodifferential operator of order —1, and this is how it is usually 
described. At least in the convex case, one may extend /3 to the boundary \rj\ = 1 of 
B*Y by fixing S*Y pointwise. Then the diagonal of B*Y x B*Y and Cbiu intersect 
in the diagonal of S*Y x S*Y, so one may view the wave front set of Fh as being 
the union of two intersecting Lagrangeans. Fh{y, y') thus has similarities to the the 
oscillatory integral operator kernels associated to two intersecting Lagrangeans of 
Melrose-Uhlmann. However, the intersection occurs at the boundary of both La- 
grangeans and is thus outside the scope of the class of Melrose-Uhlmann operators, 
presumably explaining the unusual composition law described in |Z1| . Due to the 
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explicit nature of our problem, we carry out the analysis by hand without making 

use of general operator theories. The main point to observe is that the Fourier 
integral part is of order while the pseudodifferential part is of order —1, so the 
Fourier integral part dominates. 

Proposition 3.1. Assume that dQ is smooth. Let U be any neighbourhood of 
^S'dQ- Then there is a decomposition of Fx as 

where Fi is a Fourier Integral operator of order zero associated with the canonical 
relation Cbm = graph{(3), F2 is a pseudodifferential operator of order —1 and F3 
has operator wavefront set contained in U . 

This implies an Egorov type result for the operator F\: 

Proposition 3.2. Let Ah = Op(aft) be a zeroth order operator whose symbol 
a{y,ri,0) at h = is supported away from \r]\ = 1. Put h = and let 7 be as 
above, and let (3 denote the billiard ball map on B*Y. Then 

F^AhFx = Ah + Sh, 

where Ah is a zeroth order pseudodifferential operator and \\Sh\\L^^L^ < Ch. The 
symbol of Ah is 



(3.5) 



0, B*Y. 



This is a rigorous version of the statement that Fh quantizes the billiard ball 
map. 

4. Boundary quantum ergodicity 

Quantum ergodicity is concerned with quantizations of classically ergodic Hamil- 
tonian systems. It is essentially a convexity result relating the time average 

{A) := hm ^ r U,AU;dt, 

and the (constant) space average, 

LJ{A) I 

of an observable. Here, we use the notation of the interior problem where Ut is the 
wave group. Also, u){A) is an invariant state on the algebra of observables, which 
arises from the local Weyl law: 



(4.1) lim -i- ^ {Auj,Uj) = uj{A). 

A^oo N{X) 

The system is quantum ergodic if 



(4.2) {A) = oj{A) I + K, with J^WUxKUxWhs ^ 0, 
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where || • \ \hs is the Hilbert-Schmidt norm and where Il\ is the spectral projection 
onto the span of eigenfunctions of \/A of eigenvalue < A. In terms of the L^- 
normalized eigenfunctions Uj, QE says: 

^ \{Au,,u,)-uiA)f^O, 

for any observable A. Following the work of Schnirelman, Colin de Verdiere, 
Zelditch in the boundaryless case, it was proved by Gerard-Leichtnman jOLj (C^'^ 
domains with Dirichlet boundary conditions) and Zelditch-Zworski [ZZw| (general 
case with corners) that domains with ergodic billiards are quantum ergodic. 

4.1. Proof of Theorem [ij The main steps in the proof are: 

• The local Weyl law; 

• Analysis of the classical limit state and the ergodic theorem on the 
classical level. 

• A convexity inequality to convert quantum ergodicity to the classical 
ergodic theorem. 

4.1.1. Local Weyl law. It states: 

Proposition 4.1. Let Ah he a semiclassical operator of order zero on dfl. Then 
for any of the above boundary conditions B, we have: 

(4.3) ^lim J2 i^h,u),u)) = ub{A), h, = Xj\ 

where ljb is the state defined in the table in Section 

In the case of multiplication operators, the Weyl law was first proved by S. Ozawa 

m 

4.1.2. Mean ergodic theorem for the classical limit states. It is obvious that the 
states 

(4.4) p){A) := {Au^u],u'^ 
are invariant for Fx. : 

(4.5) p'^iF^^AFx,) ^ p',{A) 

Here, as above, hj — XJ^ . 

Since any average or limit of averages of these states will be invariant, the local 
Weyl law implies: 

Corollary 4.2. The states lob are invariant under F\: ujb{FxAFx) ~ lub{A). 

The formulae for the limit states in the table in Section |0] are found by explicit 
calculation of traces, e.g. of the dual heat or wave traces. The limit measures may 
be understood as follows: 

• The Neumann limit measure is the projection to _B*J7 of the interior 
Liouville measure restricted to the set S^d^ of inward pointing unit vectors 
to 17 along dn under the map Tr{x,^) = (a;,^"^) taking a (co)-vector to its 
tangential component. 



10 



STEVE ZELDITCH 



• The Dirichlet limit measure j{q)da arises because boundary values of eigen- 
functions involve the normal derivative in both factors before restricting to 
the boundary. The symbol of hd^^ restricted to the spherical normal bun- 
dle, and then projected to B*dil is equal to 7, so we expect to get the 
square of this factor in the Dirichlet case compared to the Neumann case. 

We now consider the mean ergodic theorem for these classical limit states. Let us 
first recall the result for the standard 'Koopman' operator associated to the billiard 
map: 

T : L^B*dn, da) -> L^B*dn, da), T/(C) = /(/3(C)). 

From the invariance it follows that T is a unitary operator. When /3 is ergodic, 
the unique invariant L^-normalized eigenfunction is a constant c, and one has the 
mean ergodic theorem 

(4.6) hm ||lf;r;?(/)-(/,c)||^0. 

n=l 

But the transformation provided by Proposition l3.2l f Egorovs theorem) is defined 
by: 

This T is not unitary on L'^{B*dVt,da). The invariance ujb{A) = LUB{FlAFh) 
implies that the associated measure d^s is invariant under 

7(C) 

Simple calculations show: 

• (i) The unique positive T*-invariant density is given by 7~^d/i. The unique 
positive T-invariant density is given by ^d^. 

• (ii) T is unitary relative to the inner product ((, )) on B*dil defined by the 
measure dv = j~'^d^. 

When /3 is ergodic, the orthogonal projection P onto the T-invariant L^-eigenvectors 
has the form 

P{.f) = {{frd) = [f fii-^dy^Bh 

JB'dn 

Thus, P{aA) = [Ji^,g^aAl~^diiBh = ^^5(^)7. 
The mean ergodic theorem thus says: 

JV-l 

-Y^T'^aA^u^siA)^. 
4.1.3. Convexity. To show that 

{Au),y!'^)^ujB{A), 
along a density one subsequence of integers j is essentially to show that 

(4.7) linisup^ J2 \{iA-u;BiA))u),u) 
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Due to the novel form of the local Weyl law and the Egorov theorem, we first 
prove an auxiliary result of this kind for the quantization of the invariant symbol 



(4.8) 



with c a normalizing constant. 

Lemma 4.3. For all e > 0, there exists a pseudodifferential operator R of the form 
i4-S}) such that 

1 I 2 

A^oo JV^AJ ^^^^1 

We prove this intermediate step by the usual time-average + convexity: 



= Ea,<a {{{Ar.)N-u;BiA)R}u),u)) 
< CEx,<x{\{Ah)N^u;BiA) R\'u>;,u'^), 



where (A/Jat = ;^ Ef=i ((^/i )*^^^-^^ By the local Weyl law, the limit equals 
lubHAn -iUBiA)R)^), or 



(4.9) 



B'dn 



{crA)N - i^B{A)a{R)) dfiB 



By Proposition the symbol cr^^^^ of {A) at is J2k=o "^f^^a)- By the mean 
ergodic theorem, this converges to 

P{f^A) = X / ^Al^^ da = cujb{A) 7(g) . 

This is approximately equal to the symbol of lob{A)R. Thus, 



B-dn 



ctan - ^B{A)cr{R) \ dfiB 



becomes small as N ^ oo. Thus, j^jj^ J2\j<\ {(A — iLiBiA) R)u^j, u' 
small as iV ^ oo, proving Lemma 14.31 
Finally, we need to go from 

^ J2 \{{A~u;b{A)R^,u\ " 



is arbitarily 



iV(A) 



Xj<X 



to 



^ E \{{A-^b{A)I)u],u]) 



Xj<X 



It suffices to show that -j^j^ J2xj<x ii^ ^ ^) 
this is the case A — I above, since lob{I) ~ 1- 



is arbitrarily small. But 
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4.1.4. Non-convex domains. We now adapt the argument to nonconvex domains. 

The problem with the argument above is that the canonical relation of F\ is larger 
than the billiard relation, since it relates points on the boundary which are con- 
nected by a straight line even if the line passes outside the domain (ghost orbits). 
However, one can modify Fx so that it leaves the boundary traces invariant and so 
that its wave front relation has an arbitrarily small measure of ghost orbits. That 
this should be possible was emphasized to us by M. Zworski. 

As is easy to see, the property F\.v!j = (which follows from Green's formula) 
is valid for any choice of Green's function on a neighborhood of il. So we modify the 
metric outside the domain, while keeping it Euclidean inside. Let 6 be a smooth, 
compactly supported nonnegative function on M" which vanishes on fl. Consider 
the metric 

gs = {l + s6)5fEuclidean On M". 

For sufficiently small s, no geodesies of g,, starting at a point in have conjugate 
points in Let Gs{z,z,X) = (Ag — (A + iQi)'^)~^{z, z') denote the kernel of the 
outgoing resolvent of the Laplacian on R" with respect to the metric Qs- It has the 
parametrix 

Gt{z,z',\) = A"-2e*^'^^"*^(^'^')a(^,^',A), zj^z'. 

We then define: 

F',{y,y')^2^Gsiy,y',X), y,y' eY, 
and average over s to obtain the operator 

(4.10) Fx= j x{s)Flds. 

Jo 

Here x is smooth, supported in (0, 6), and nonnegative with integral 1. This averag- 
ing removes all but an arbitrarily small measure of ghost orbits or spurioius points 
from the WF'{Fx), i.e. pomts{y,r],y' ,r]'), with t] = dy\y - y'\, tj' = -dy>\y - y'\ 
such that the line yy' leaves O. We still have 

and in addtion the averaged operator has, as in the convex case, a decomposition 

where Fi is a Fourier Integral operator of order zero associated with the canonical 
relation Cbui, F2 is a pseudodifferential operator of order —1 and F3 has operator 
wavefront set contained in U. Moreover, the symbol of A, a is the same as for 
Fx- Thus, the proof quantum ergodicity in the non-convex case now runs as in the 
convex case, with this averaged operator in place of Fx- 

5. Norm estimates of boundary traces 

We now turn to the results on norms of boundary traces. They are derived 
from an analysis of the singularities of the boundary trace 

00 

(5.1) i;|(t,g,g)=^e^*^^>5(g)|2 

j=i 

of the kernel of the wave operator cosf-v/Ag on O with boundary condition B. 
The nice feature of the boundary trace of EB{t,x,x) is that the singularity at 



BILLIARDS AND BOUNDARY TRACES OF EIGENFUNCTIONS 



13 



t = becomes uniformly isolated from other singularities, while the interior kernel 
E{t,x,x) has singularities at t = 2d{x) arbitrarily close to i = 0. We note that 
E%{t, q, q') is a spectral transform of the Dirichlet-to-Neumann kernel. 

We recall that the boundary traces are normaHzed by 11^^11^2(0), so that it 
is of interest even to obtain estimates on ||w^||L2(ao)- In H2.3(l - (|2.4|l we obtained 
asymptotics of a density one sequence of boundary traces in the ergodic case. In 
general, the results are: 

• Under a non-trapping assumption, CiA < ||u^||L2(ao)/||i*j||_L2(Q) < C2A in 
the Dirichlet case IBlHIIHT] : 

• For any smooth domain and metric, ||uJ||L2(aQ)/||wj||^2(f2) = 0(A^^^) in 
the Neumann case (Tataru, T , Theorem 3); 

• For a Sinai billiard on a torus (the exterior of a convex smooth domain 
in the torus), ||uJ||L2(aQ/||wj||i2(o) — 0{\^J^) in the Neumann case 
Theorem 5). 

We introduce some notation. Given (g,?]) G B*d^, we let ^{q,ri) e >5'*„ be 
the inward pointing unit vector to Q obtained from {q, rf) by adding a multiple of 
the unit normal. We further denote by $* the broken bicharacteristic flow of the 
wave group, i.e. the flow which carries singularities of the solution of the wave 
equation. Finally, we let denote the tangential projection of t; G S*Vl at q e dO,. 
Finally, we denote by 7^ the boundary trace in the q variable corresponding to the 
boundary condition B. 

5.1. Singularity of boundary trace of wave kernel at t = {]. Isolation of the 
singularity at t = of Eg{t, q, q) follows from simple wave front set considerations. 
We first note that 

WF{j^j§E{t,q,q')) C {{t,T,q,rj,q\7^') : ^(9, r;))]^ = {q',rj'), r = 

which follows from propagation of singularities for the wave equation. Composition 
with the boundary trace just pulls back (i.e. restricts) this wave front relation to 
the boundary. 
It follows that 

WF{^!^^!^,E{t,q,q)) C {(i, r, g, ,7, g, 77') : [<f'{qA{q,v))f = {q,v'), r = -|e(g,??)|}. 

Thus, the singularities of the boundary trace lqlqiE{t, g, g) at g € dO. to broken 
bicharacteristic loops based at g in f2. Let us first consider loops in convex domains. 
They are either: 

• closed geodesies onY — 9ri; 

• m-link transversal reflecting rays with m vertices in Y , one of which is at 
g and which satisfy Snell's law of equal angles except possibly at g. 

In the non-convex case, the description is more complicated since there exist 
additional gliding rays and since boundary geodesies need not carry singularities. 
For instance, for Euclidean plane domains in dimension 2, bicharacteristics enter 
and exit the boundary at inflection points and glide only over the convex part of 
the boundary. 
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By definition of normal singularity, there exists r e such that t^E''{t, q, q) e 
C°° near < = 0. Thus, there exist coefficients aj{q) such that 

oo 

(5.2) i?^(i,g,g)^t-'-^a,(g)t^ 

3=0 

The next step is to calculate the coefficients of the singularity. 
Using a method of Ivrii, we first prove: 

Proposition 5.1. The singularity of E''{t,x,x) at t = is a classical conormal 
singularity. 

Granted the proposition, we can calculate the coefficients (and the order) of the 
singularity using Hadamard type variational formulae for eigenvalues: 

5\j^ j^p\u){y)\-'dA{yl 

where 5 denotes a variation of either the boundary or the boundary conditions and p 
represents the tangent vector to the variation. Thus, we can determine the integrals 
p(q)aj{q)dA{q) and hence the coefficients aj{q) asymptotically by considering 
the variation of the wave trace 

3 3 3 ^ 

5.1.1. Dirichlet. In this case, we vary the boundary in the normal direction with 
variation vector field pv. The wave trace formula has the form: 

^e'*^^ = CnVol,,{fl){t+iO)-"+Cn-iVoln^i{dn){t+iOy+'^+ai{t+iO)-"+^+--- 

3 

where • • • represents lower order terms (cf. I , Theorem 2.1 or |^, Corollary (3.6)). 

It follows that 

(5.4) 

Ee'*^' J p|u5(y)|'dA(y) =C„5Fo?„(r!)(t+zO)-"^i+C7„-i<5Fo/„_i(51))(t+iO)-"+--- 

Here, we divided by the coefficient it in H5.3I) . The variation of Vol{il,) is non-zero, 
so we get a {t + iO)~"~^ term on the right. 

5.1.2. Neumann boundary conditions. In this case, we vary the boundary condi- 
tions, and therefore consider more general boundary conditions of the form 

d,cj){q) + n{q)4>{q) = 0, g G d^. 

We denote by 5 the first variation relative to a change in the boundary condition 
from K ^ K -\- ep. By the Hadamard variational formula (studied by S. Ozawa for 
general boundary conditions), the jth Neumann eigenvalue pj has the variation 

5p, = f \u''^{q)\'p{q)da{q). 

We now consider the first variation of the Neumann wave trace under variations 
of pseudodifferential boundary conditions. We obtain: 

(5.5) 6TrE{t) ^itY.Spje''"^. 

3 
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The p term only influences the third term of the singularity trace expansion at t = 
for Neumann-Robin boundary conditions (see |GM| for the calculation) Hence the 
variation has the form 

(5.6) 5TrE{t) = C„( / pdp){t + ^0)""+^ + • • • , 

J S'dO. 

from which we conclude 

(5.7) Ve**^^[/ p\u]\HA\=Cn{l pdfi){t + 10)-"+' + ■ ■ ■ 

^ Jan Js'on 

5.2. Spectral asymptotics of boundary traces. We apply standard Tauberian 
theorems to obtain spectral asymptotics of boundary traces from the singularity of 
the boundary trace of the wave kernel. 

Proposition 5.2. We have: 

( CA"+2 +0(A"+i), Dirichlet 

j:x,<\ [ CA" +0(A"-i), Neumann. 

Proof. It follows by Proposition l5 . II and a standard Tauberian argument that there 
exists a two-term asymptotic expansion for some principal coefficient. We may 
determine it from an integrated version of the asymptotics where we integrate the 
left side against a smooth function p on dfl. The integrated version follows from 
()5.4|l in the Dirichlet case and (|5.7|) in the Neumann case. 

Lemma 5.3. For any smooth p on dfl, 

I- [ /g^pdAA"+2 +0(A"+i), Dirichlet 

E / p\u]{q)\'dA={ 
r.\j<\''^ y ^Q^pdAy + 0{X'-'), Neumann. 

□ 

We now improve the result in the generic case: 
Theorem 5.4. Suppose that the set of loops at q has measure in B*dfl. Then 

( CA"+2 +o(A"+i), Dirichlet 

E W\' = \ 

j:Xj<\ [ CA" + o(A"' ^), Neumann. 

It follows that I l^oo^gQ) ~ o(A~5~) in the Dirichlet case, resp. o(A~2— ) in the 
Neumann case. These results are sharp. 



Proof. The proof follows the outline of that in |SZj . We write the left side as the 
boundary trace E^^ ^-^ {q, q) of the spectral projections, and define the remainders in 
the local Weyl laws by 



^[0,A] (?> Q) 



CA"+2 + i?c(A, g), Dirichlet 
CA" -I- i?Ar(A, q), Neumann. 
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We first show that if the set of billiard loops a.t q E dfl has measure zero in 
*dn, the 
for A > A, 



B*dil, then given e > 0, we can find a ball B centered at q and a A < oo so that 



(5.8) 



\RD{X,q)\ <eA"+i, qeB, 



|i?jv(A,g)| <eA"-i, qeB. 

To prove this, we need some more notation. The loop-length function on B*dfl 
is the lower semi-continuous funciton defined by 

{mi{n > : no f3"{q, -q) = q}, 
oo, if no such n exists. 
The set of loop directions at g G dO, is defined by: 
(5.10) Cq = {ij^Bldni/L*{q,ij)^Q}. 

We construct two semiclassical pseudodifferential operators b{q, D), B{q, D) ~ 
I — b{x, D) on L^{dU,) with the property that B is microsupported in the set where 

[0,A]( 



L*{q, Tj) » T. We then study Et ^, (q, q) by writing it in the form 



(5.11) ^o,,] [q, q, )^[{B + h)Ei{B + b)] [q, q). 

We choose b e C°°(S'"-i) so that 



(5.12) / 6(77)da(77) < l/T^ 
and 

\L*{q,r])\ <1/T, onTVxsuppB, 
where A/" is a neighborhood of q. We then put 

5(0 = 1-6(0. 

By construction, 

(5.13) E\t)B*{q,q), BE''{t,q,q)eC^{Q,T). 

A calculation using the conormal sing ularity at t = of BE''{t, q, q), bE^{t, q, q) 
shows that, for A > 1, 

(5.14) |i?Jv(A,q)| < CT-^X"--^ +CtX"--'^, 

where Ct depends on T but C does not. This of course yields ()5.8|l . The argument 
is similar for the Dirichlet case. 

□ 

5.3. Applications to eigenfunctions. We now use that bounds on eigenfunctions 
(or for spectral projections for intervals of shrinking width) can be obtained from 
the jump in the remainder: 

(5.15) W]{q)\^^R{\,q)-R{\~Q,q). 

j:\j=\ 

To complete the proof, we observe that for fixed q £ dfl and any s > then one 
can find a neighborhood Afe{q) of q and an A£(g) so that when A > Ai^^xq) and 
y e Afeiq) we have \R{X,y)\ < eA""^ This implies that \u''^{y)\ < eA^""!)/^ if 
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y G Me{q) and A > h.e{<l)- Since M is compact and since the Af^{q) form open cover 
of M, we may choose a finite subcover and extract the largest A^{q). For this A^, 
we get 

K(2/)| < £A("-')/^ A>A„ 
The o(A("-i)/2) bound follows since depends only on e. 

5.4. Examples. For generic metrics and/or boundaries, there are no recurrent 
points jSZ| . Moreover, convex analytic domains in R" never have have recurrent 
points on the boundary. Suppose to the contrary that there exists xq G such 
that a positive measure of geodesic rays starting at xq return to xq at the same 
time. By analyticity, all rays starting at xq return to xq, and they must all return 
at the same time. In particular, boundary geodesies have to return to xq. So do 
creeping rays which make a small angle to the boundary, and which accumulate at 
boundary geodesies. But the creeping rays which approach a boundary geodesic 
are strictly shorter than the limiting boundary geodesic in the Euclidean metric, so 
the return time could not be the same. Therefore, the L°° estimate on boundary 
traces of eigenfunctions or spectral projections is never sharp for convex analytic 
Euclidean domains. 

On the other hand, the L°° bound is sharp for the Euclidean half-circle or for any 
sector of a circle. Indeed, the invariant Neumann eigenfunctions (under rotations) of 
the disc achieve the interior L°° bound at the center. On a half-circle or sector, they 
remain Neumann eigenfunctions and their boundary traces achieve the maximal 
L°° bound above. In the case of Dirichlet boundary conditions, the bounds are also 
saturated on such domains (by taking the boundary trace of the imaginary part of 
the eigenfunction transforming by e'^ under rotation by angle 0). 

Above, we have mainly considered simply connected domains, but the same 
questions may be posed for multiply connected plane domains (for instance). It 
seems likely that any analytic (or even piecewise analytic) domain which achieves 
the maximum L°° bound must be simply connected. This would generalize the 
topological result of |SZj that the maximum L°° bound can only be achieved for 
metrics on the sphere. 

The question thus arises to find the maximal growth rate of L°° norms of bound- 
ary traces of eigenfunctions of smooth Euclidean domains, (or of convex analytic 
Euclidean domains), and which domains achieve the bounds? In fact, the answer 
depends on whether one studies relative sup-norms ||M^||L°c(an)/||Mj||L2(gQ) or ab- 
solute sup-norms ||M5l|L°°(aa)- Moreover, the sharpness depends on whether we 
consider individual eigenfunctions or spectral projections for shrinking intervals. 

We do not know the answer to the question, but offer some speculations. By 
recent results of Smith-Sogge jSSl) following earlier results of Grieser (in his un- 
published PhD thesis and in (see also [H])), it is known that whispering gallery 
Neumann modes of two-dimensional convex Euclidean domains saturate the 
norms with 2 < p < 8 because they live in thin A^^^^ layers around the boundary 
and therefore are of size X^^^ there. The analogous problem in higher dimensions 
is still open. 

However, they are not extremal for relative sup norms because are spread out 
all over the boundary. For instance, in the case of the disc, the relative sup norm 

equal to one. A better guess is that boundary traces 
of modes (or quasimodes) associated to stable elliptic orbits are extremals for the 



18 



STEVE ZELDITCH 



relative sup norm ||u^||Loo(an)/||Wj||L2(af2'). They live in A ^1"^ tubes around the 
orbits and therefore the modes have size A'"~^^/^ along the orbit in dimension n. 
The boundary trace is therefore concentrated in small balls of radius around 
the bounce points of the orbit. So the norm of the boundary trace should be 
- 1. 
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